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ABSTRACT 

We present results of a study of very low mass halo stars. Using a sample of proper 
motion stars identified from plate material taken as part of the first and second Palomar 
Sky Surveys, we measure the space density, metallicity distribution, and kinematics of 
the Population II M subdwarfs. Our overall luminosity function is in good agreement 
with previous analyses of the space density of nearby very-low-mass halo subdwarfs, 
and confirms the discrepancy between local analyses and the space densities inferred 
from deep HST starcounts. We show for the first time that both the metallicity 
distribution and kinematics of late-type halo subdwarfs are consistent with those of 
their higher mass metal-poor counterparts. Dividing our sample by abundance, we find 
no evidence that the mass function of field halo stars is dependent upon metallicity. 
We provide data for three nearby subdwarfs that may merit additional observations. 

Subject headings: Galaxy: halo — stars: low-mass, brown dwarfs — stars: luminosity 
function, mass function — stars: Population II 

1. Introduction 

Long-lived low-mass stars are a local relic of the earliest stages of the Galaxy's formation. 
Dubbed the Population II in contrast to the disk Population I, these stars make up the Galactic 
halo (sometimes called spheroid). Although comprising only a very small faction of the stars in 
the Galaxy, Population II stars provide one of the most important records of Galactic history. 

The Population II luminosity function is of particular interest since its derivation is a 
requirement in most studies of the stellar mass function. Determination of the stellar mass 
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function as function of Galactic epoch and metallicity provides significant constraints upon the 
theory of star formation. Evolution of the mass function with time or metallicity might be 
expected since conditions in star-forming clouds were likely much different in the protoGalaxy 
(c.f. Zinnecker 1995| ), Moreover, since halo brown dwarfs have faded far below the threshold for 



even near-infrared detection, the form of the subdwarf mass function near the hydrogen-burning 
limit represents the most effective technique (in the absence of gravitational lensing detections) 
for estimating the Population II dark matter contribution. 

Despite the intuitive simplicity of the process of counting stars to determine the luminosity 
function, there have substantial disagreements between different studies in both the normalization 
and slope of the luminosity function (see the review of Mould 1996| ). Two recent surveys have 



measured the space density of very-low-mass Population II stars, the M subdwarfs. Counts of 



faint, distant M subdwarfs using HST (Gould et al. 1998) predict ~ 2.5 times fewer stars at 



My = 11 than are found locally by Dahn et al. (1995). Dahn et al. (1995) base their study upon 
a sample of the local high proper motion stars in the LHS Catalog ( [Luyten 1979| ) for which they 
have measured accurate trigonometric parallaxes. 

In this paper, we determine the luminosity function of the field Population II M subdwarfs 
using a new proper motion survey based upon plates taken as part of the first and second Palomar 
Sky Surveys. In Section we discuss some general problems with the selection of halo stars. In 
Section |j| we describe the data used in our survey and the methods used to select a halo sample. 
In Section we use our data to measure the M subdwarf luminosity function. In Section || 
we compare our luminosity function to other published luminosity functions and discuss its 
significance. In Section ^, we summarize our results. In the Appendix, we describe some stars that 
may be of particular interest. 



2. Selection of Halo Stars — General Principles 

Measurement of the Population II luminosity function is considerable more difficult than the 
measurement of the Population I luminosity function, for the simple reason that the Population 
II represents a tiny fraction of stars locally. One must either count very faint stars at large 
distances from the Galactic plane (Gould et al. 1998), in which case one sacrifices the ability to 



obtain spectroscopic data, or else determine some way to identify the local Population II stars 
by excluding the local ~ 99.8% of stars that belong to the disk. The latter course requires an 
efficient method of eliminating the disk stars before time-consuming followup observations are 
obtained. As members of a high velocity-dispersion, low rotation population, subdwarfs tend to 
have high heliocentric velocities, and choosing an initial sample of proper motion stars increases 
the halo contribution from ~ 0.2% to ~ 30% (Schmidt 1975). Such selection criteria introduce 
some kinematic bias, but lead to a much more manageable sample. Nevertheless, such a sample 
is still dominated by old disk (Population I) and "thick disk" (Intermediate Population II, IPII) 
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starsf]. 

As recognized by Schmidt (1975) and further discussed by Bahcall & Casertano (1986, 
hereafter BC), even a small proportion of contamination by high velocity disk stars can lead a gross 
overestimate of the space density of the halo. In order to ensure a pure halo sample for his study, 
Schmidt imposed a restrictive velocity criterion (vtan > 250 km/s) and then applied a correction 
to account for halo stars with smaller velocities. BC used Monte Carlo simulations to estimate a 
correction factor of 3.03 for a cutoff velocity of 220 km s _1 (this correction is dependent on the 
kinematics assumed for the underlying population). Most recently, this technique has been used 
to derive a halo luminosity function from the LHS catalog stars with 0. 8"yr _1 < M < 2.50"yr~ 1 
and 11.0 < m r < 18.1 (pahn et al. 1995| ). 

If trigonometric parallaxes are available for the target stars, then the distances and tangential 
velocities are estimated easily and the appropriate kinematic criteria can be applied. This is 
the case for the Schmidt (1975) sample and the Dahn et al. (1995) samples. Measurement 
of a useful trigonometric parallax requires years of observing — and the the target must be 
relatively close (d ^ 100 pc). Expansion of the sample to cover a larger space volume requires 
a less direct method of distance estimation. Photometric parallaxes provide that method, with 
the calibrating color-absolute magnitude relations provided by subdwarfs with high quality 
trigonometric parallaxes. 

The efficiency of the initial selection of halo candidates can be increased by use of the 
"reduced proper motion", H, first used by Luyten (1939). It is defined as 

H = m + 51og/i + 5 = M + 51og|^ (1) 

where m is the apparent magnitude and M is the absolute magnitude. Note that H is distance 
independent and determined entirely by observables. If we assume a color-magnitude relationship, 
we can plot constant tangential velocity contours on the data. Figure [l] illustrates the concept, 
where data are plotted for 2111 M dwarfs within ~30 parsecs of the Sun(Reid et al. 1995; Hawley 
et al. 1996). (Note that the increasing proper motion bias for redder stars in the Gliese &: 
Jahreiss 1991 nearby star catalog is clearly evident.) The lines illustrate where a disk population 
with halo-like velocities will lie, and how halo stars will be even more easily distinguished in the 
diagram. 

As subluminous stars in the (M/, (R-I)) and (M/, (V-I)) planes, halo stars have systematically 
larger H than disk dwarfs at a given color, mimicking the effects of increasing vtan- In contrast, 
the coolest metal-poor subdwarfs are superluminous in the (My, (B-V)) plane (Gizis 1997), and 
a cutoff in the (Hy, (B-V)) diagram tends to discriminate against lower-abundance stars. We 



3 The term Intermediate Population II, used instead of "thick disk" or "extended disk," should not be confused 
with the Population II halo. Considerable evidence exists that the two populations are kinematically and chemically 
discrete (Majewski 1993). 
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have therefore defined our sample of candidate halo subdwarfs as those stars with reduced proper 
motion, Hj, which exceed the value predicted for disk dwarfs with vtan = 220km s -1 and of the 
appropriate (V-I) or (R-I) color. Since the halo stars are subluminous, this selection corresponds 
to a cutoff velocity that is a function of metallicity - an effect which is taken into account in the 
subsequent analysis. 

3. Selection of Halo Stars — Data and Method 

We utilized the techniques described above to define our sample. Proper motion stars were 
identified using the first and second epoch 48-inch Schmidt photographic plates, as described 
in (Section |3.l[) . An initial cut of the sample was made to select large tangential velocity stars 
(Section |3.2| ). We then obtained spectra of each candidate halo star to estimate its metallicity and 
radial velocity (Section ^ ). 

3.1. Proper Motions 

Table |l| lists the fields analyzed in the current survey. In each, the effective area is ~ 25 square 
degrees, with most having plate material spanning a baseline of at least 40 years. In each case, 
the first epoch plates were taken with the 48-in. Oschin Telescope as part of the first Palomar 
Observatory Sky Survey (hereafter POSSI). For the most part, the second epoch plates were taken 
in the 1980s by either the 48-in. Oschin Telescope (in the course of the Palomar Observatory 
Second Sky Survey, POSSII) or by the U.K. Schmidt Telescope (UKST). These second epoch 
plates were used by Tinney, Reid, & Mould (1993, hereafter TRM) in their determination of 
the disk luminosity function. All of the POSSI plates were scanned by the APM facility at the 
Institute of Astronomy, Cambridge, while the POSSII/UKST plates were scanned by COSMOS 
at the Royal Observatory, Edinburgh. One additional field, centered near the North Galactic 
Pole, was included. As described by Reid (1990), the second epoch plates for this field were taken 
at Palomar in 1976, and all plate material was scanned by COSMOS. Each of these scans yields 
positions (both pixel x,y positions on the plate and right ascension, declination positions on the 
sky), magnitudes, and morphological classification as stars, galaxies, "merged stars" or noise. All 
objects classified as noise were excluded from further analysis. 

As a preliminary to identifying proper motion stars, all objects (stars and galaxies) identified 
in the scans of the plates taken at the same epoch must be matched. In the case of the first 
epoch POSSI data, this is straightforward, since the O and E plates (corresponding roughly to 
photometric B and R observations) were taken on the same night. Thus there is no real motion 
between the time of the observations and the plate scans can be matched by demanding positional 
coincidence to a high level of precision. For our second epoch TRM data, the IIIaF (photometric 
R) and IV-N (photometric I) plates were not observed on the same night. TRM paired the objects 
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identified on their scans using a 3 arcsec box. Since the TRM plates were taken ~ 3 years apart, 
this procedure loses the highest (/x ;> 1 arcsec yr -1 ) proper motion stars. In our survey, this is 
not a concern because the 40 year separation between the first and second epoch, with no plates 
at intermediate epochs, prevents our identifying reliably stars with high proper motion due to the 
large number of possible pairings. The effective upper limit for our survey is /i=0.375 arcsec yr _1 , 
and that limit is taken into account in our luminosity function analysis. 

We should note that the TRM detections are restricted to the region within 3 degrees of 
the plate center in order to avoid the vignetted regions near the edge of the plates, and that the 
COSMOS scans did not cover the entire plate. The appropriate area of the survey is given in 
Table |]. Moreover, the the POSSI field centers do not correspond to the POSSII/UKST field 
centers. Thus, between one and four POSSI plates are paired with the TRM fields, as listed in 
Table @. 

The identification of the proper motion stars was made as follows. First, an initial pairing of 
objects within a 3 arcsecond box (A < 3") was made. In each field, the surface density of celestial 
objects is sufficiently low that matching within this radius results in an unambiguous pairing. 
Figure ^ illustrates the agreement between the rp and E magnitudes for one representative pairing 
(POSSII Field 513 and POSSI E102). Most of the outliers evident in this plot are correct pairings 
of extended objects (galaxies or merged stars) — they are outliers due to differences in the 
definitions of magnitudes used by the two machines, as well as the different plate characteristics. 
For this study, we are interested only in stars, so such outliers are unimportant. Having paired the 
unambiguous matches, we obtain a list of unmatched objects — these objects have no counterparts 
within 3 arcseconds at the second epoch and are good candidates for proper motion stars. 
Eliminating paired objects, we find all possible matches within 15 arcsec of unmatched objects. 
We exclude pairings which have magnitudes inconsistent with the E -rp relation delineated by 
the close matches — inconsistent being defined as a discrepancy of from >1 magnitude rp = 12 to 
> 1.5 magnitudes at rp = 19. These limits are "conservative" in the sense that some false pairings 
are included in order to avoid excluding true proper motion stars. 

The vast majority of objects on the plate should show insignificant proper motion between 
the two epochs, but directly matching the right ascension and declination catalog positions leads 
to many (false) proper motion objects. Those spurious motions arise from systematic errors 
in the relative positions of the two catalogs (caused by differing plate centers, inadequacy of 
the plate model, plate variations, telescope differences, etc.) These systematic positional errors 
can be corrected by determining the transformation between the original pixel positions for the 
two epochs (a?i,yi and £2,2/2)- Both our own experience with the TRM fields and the Reid 
(1990) analysis showed that a polynomial solution over the entire plate works poorly, producing 
many false proper motion stars. Analysis over smaller regions, however, is effective. Each plate 
was divided into approximately one by one degree regions and the positional transformation 
determined using second-order two-dimensional polynomials. Only stars with A < 3" matches 
were used for the initial co-ordinate transformation, iterating twice to exclude objects with large 
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(> 2") residuals. We then used these polynomial solutions to determine the proper motions of all 
stars in the field. The stars that prove to have /j, > 0.1 arcsec yr" 1 are included in our sample of 



proper motion stars. The selection of the candidate halo stars is described in Section 3.2 



These proper motions are measured relative to faint stars, which have typical distances 
of ~ 500 — 4000 pc. Thus we expect the reference frame to have a small, but non-zero, mean 
proper motion. We derive the conversion from this relative frame to an absolute (extragalactic) 
frame using measurements of objects classified as galaxies. The latter are identified using the 
^-parameter defined by Picard (1991). We take 200 < (p < 600 for these purposes, which provides 
a sample dominated by galaxies. The upper cutoff is taken to exclude the "fuzziest" galaxies - 
we found that these "fuzzy" galaxies have large spurious proper motions, which we attribute to a 
combination of the different plate material used at the different epochs, and the different scanning 
and centroiding techniques used at the APM and COSMOS. In Figure |3[ we plot the distribution 
of galaxy proper motions and stellar proper motions for the fields. It is clear that the galaxies do 
indeed show a systematic shift with respect to the stars. Table || lists the median proper motions 
(A/i) derived for the galaxies, which at most amount to 5.2 mas yr _1 . We have checked these 
values by comparison with the projection of the solar motion and rotational lag for the IPII and 
halo populations. The direction of motion derived by the galaxies (i.e., the sign of the proper 
motions) agrees with the expectations. These values given in Table | are used to correct the 
relative proper motions to an absolute system: 

Habs = l^rel ~ A// 

We also give the observed standard deviation of the galaxian motions, with values in the range 
10.6 — 17.0 mas yr _1 . These provide an estimate of the uncertainties of our proper motions, but 
we expect that the extended galaxies should have poorer positions, and hence less accurate proper 
motions, than is the case for stars. 

It is evident in Table ^ that Field 831 shows the largest random errors and has absolute proper 
motions correction that are discrepant at the ~ 2 mas level from the adjoining fields. R. Mendez 
(private communication) has kindly computed model proper motions for this direction based upon 
the galactic structure model of Mendez & van Altena (1996). The predicted mean stellar proper 
motions of (fj, a ,fj,g) = (3.8, —4.2) mas yr -1 are in excellent agreement with our observations in 
Fields 829 and 832. We take this as confirmation that our simple correction technique using 
galaxies is adequate for our purposes. The data for Field 831 are evidently of lower astrometric 
quality, but the accuracy remains adequate for identifying the high proper motion stars analyzed 
in this study. 



The NLTT catalog ( [Luyten 1979-1980 ) offers an opportunity to check both the completeness 



and accuracy of our measurements. As noted by many authors, identifying the NLTT stars can 
be difficult since some of the positions are quite poor (for example, we find a star that matches 
the magnitudes and proper motion of LP 621-64, but is 10 arcminutes east of Luyten's position). 
Excluding stars which are too bright, too faint, or too blue (i.e., too faint to be seen on the POSSII 
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IV-N plates), we find that 122 of 137 NLTT stars are recovered, indicating a completeness of at 
least 89%. We find that the standard deviations of Ha° SS ~ \^ LTT and ^f oss - ^ LTT are in 
the range 10 to 25 mas yr _1 . Luyten estimated his motions to be accurate to ±15 — 25 mas yr -1 
The combination of our uncertainties estimated from the galaxies and the NLTT uncertainties 
account for the observed scatter. 



3.2. Reduced Proper Motions 

Selection of halo candidates using reduced proper motion requires photometry in addition to 
the proper motions. We have used the photometry from the second epoch plates which provides 
us with R and I measurements for each star (except in the NGP field, for which we have V and I). 

We adopt the photometric calibration determined by TRM using CCD measurements for each 
plate. Because their ("Palomar") calibration is on a unique system, we describe it briefly here. 
A complete description is given by TRM. The best defined photometric system in R and I for 
late- type stars is the Cousins system ( Cousins 1973 )^1 However, M dwarfs have strong color terms 



in this system because the Rc* band filter extends far into the I band — as a result, the effective 
wavelength of the Rc filter moves redward for cooler stars, as discussed by Bessell (1986). The 
POSSII IIIaF emulsion plus filter combination has a sharper long-wavelength cutoff at A6900A, 
so the color term differs significantly from the standard Cousins system. The POSSII response is 
close to the Gunn r filter defined by Thuan & Gunn (1976). Therefore, TRM's CCD calibrations 
of the fields used Gunn r^* and Iq filters but observed Cousins late-type standards. These CCD 
magnitudes were matched to the COSMOS photographic magnitudes. They found that no 
significant I color term was necessary to calibrate the observations but, as expected, a strong 
R color term was necessary. Thus, their photometric observations consist of an R magnitude, 
denoted rp, an I magnitude, denoted ip, and a color, denoted (r — i)p, for each star. The relation 
between the standard Cousins color and the Palomar color is 

(r-i) P = 0.162 + 0.439(12 - I) c + 0.671(12 - lf c - 0.36(1? - lf c + 0.074(1? - lf c (2) 

ip = Ic (3) 

Note that TRM use the relation ip = Iq — 0.046, but examination of their Figure 2 indicates that 
there is no offset between ip and Ic for (R — I)c < 1.2, which is the color of the halo stars in our 
survey. The V and I photometric calibration for the NGP field is described in Reid (1990). 

The reduced proper motion diagram for Field 513 is shown in Figure |j. The solid line marks 



4 As in TRM, we use the "R" and "I" to denote generic R and I photometry on any system and indicate particular 
photometric systems with subscripts. Thus R-Ic is Cousins R-I, (r — i)p is TRM's Palomar system, and (r — i)c is 
Gunn R-I. 
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the disk main sequence denned by TRM, 

Mi = 5.56037 + 0.458615 x(r-i) P + 1.32084 x (r - if P (4) 

offset to match vtan = 220 km s" 1 . Stars falling below this line (larger H/) are identified as halo 
candidates for spectroscopic followup. Recently, we have pointed out that the disk main sequence 
shows evidence for a "kink" at spectral type M4.0 V flGizis fe Reid 1996| ; Reid fe Gizis 1997|) . 
While this feature is not included in the TRM calibration, it occurs at (R — I)c ~ 1-6, so it is not 
relevant to the current selection of candidate halo subdwarfs. 

While our tangential velocity limit is set at 220 km s _1 for solar abundance stars, M 
subdwarfs are selected at smaller v tan - The esdM enter the candidate pool for v tan > 75 km s _1 
and the brighter sdM at vtan > 125 km s~ 1 . The most likely contaminants, modestly metal-poor 
([m/H] ~ —0.6) Intermediate Population II M dwarfs, lie ~ 0.5 magnitudes below the disk main 



sequence (Gizis 1997, Figure 7), and therefore enter at vt an > 175 km s 



3.3. Spectroscopy 

We obtained optical spectra of the candidate halo stars in order to determine radial velocities 
and metallicities. The Hale 200-in. and the Las Campanas Du Pont 100-in. telescopes were used 
to obtain data in the wavelength range 6100 — 7300^4 with ~ 3A resolution. These parameters 
were chosen to correspond to the observations of nearby M dwarfs ( Reid et al. 1995| ; |Hawley et al 



19961 ) an d M subdwarfs flGizis 1997Q that we have already published. For the 200-in. observations, 



we used the Double Spectrograph (Oke and Gunn 1982). In August 1995, the blue camera was set 
to observe 6000 — 6900 A and the red camera was set to 6700 — 8000 A using 600 1/mm gratings 
blazed at 4000A and 10000A respectively. In October 1995, a new red camera was installed in the 
double spectrograph and was used in all subsequent runs to observe the region A6000 — 7400A at 
1.4 A pix" 1 with the 600 1/mm grating blazed at 10000A. With this setup, the blue camera was 
used to cover 4000 — 5500A with a 300 1/mm grating blazed at 3990A In practice, the latter data 
were useful only for the bluest ((r — i)p < 0.6) stars. The Du Pont observations used the modular 
spectrograph with a 1200 line grating blazed at 75004. 

From these spectra, we determined radial velocities and spectral types. During each observing 
run, we observed K and M radial velocity standard stars drawn from the list of Marcy & Benitz 
(1989). A velocity for each candidate halo star was determined by cross-correlating with the best 
match standard star. Comparison of subdwarfs with known radial velocities observed during the 



same runs shows that the accuracy is ±20 km s _1 ( Gizis 1997 ). Using these radial velocities, we 
measured TiO and CaH bandstrength indices, as defined originally by Reid et al. (1995). Each 
index measures the ratio of the flux (F u ) within the absorption feature to nearby pseudo-continuum 
points. Table || lists the wavelength definitions of the features (W) and the pseudo-continuum 
(S1,S2) points. We have defined a spectral classification system using these indices ( Gizis 1997[ ). 



Stars are classified as disk stars (M V), M subdwarfs (sdM), or extreme M subdwarfs (esdM). We 
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have shown on the basis of model atmosphere calculations ( [Allard &: Hauschildt 1995| ) that these 
empirical classifications correspond to [m/H] ~ —0.5, [m/H] ~ — 1.2±0.3, and [m/H] ~ — 2.0=L0.5 
(where [m/H] ~ [Fe/H]; Gizis fc Reid 1997| ). Each halo candidate has been classified using this 
system. 



4. The Luminosity Function 

In this section, we calculate the I-band luminosity function from the data described in 
Section ||. The luminosity function, $(Af/), is defined as the number of stars per cubic parsec 
per magnitude in a bin centered at Mj . We use Schmidt's (1968) 1/V max technique to estimate 

$ obs = Y,T7— (5) 

Vmax 

Vmax = -^(d max — d^iri) (6) 

Where O is the solid angle on the sky of the survey, d max is the maximum distance that a given 
star could have been detected in this survey, and V ma x is the corresponding volume. A minimum 
distance, D m i n , also appears to account for the distance at which the star would have evaded 
detection due to the upper proper motion limit of fi ma x = 0.375" yr" 1 . The maximum distance 
may be limited by either the lower proper motion limit (n m in = 0.100" yr -1 ) or the limiting 
magnitudes (Rq,Io), hence for a star, distance d, magnitudes (R, I), 

d max = d x Min dex[0.2(R - R)]; dex[0.2(I - I)]) (7) 

V /^min / 

We estimate the associated error by assuming that each star contributes an uncertainty of l/V max 
QFefton 1976Q . Then 



4 = E^- ( § ) 

' max 

As noted in Section ^, we apply a tangential velocity cutoff in order to ensure a "pure" halo 
sample. We adopt a value of v cu t ff = 200 km s _1 in our main analysis. As a result, a correction 
must be applied to the apparent luminosity function in order to obtain the space of density of halo 
stars. Thus, 

<S>halo = -$obs (9) 

X 

The discovery fraction x is a function of position on the sky and the adopted v cu t ff- In our case, 



we compute x f° r each field in the survey. This procedure is described in Section \L2 
A check on the completeness of the survey can be made by considering 

V \ If d x3 

Vrr 




(10) 
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which considers the ratio of the volume V corresponding to the star's actual distance d to V max . 
In the case of a spatially uniform sample, such as the halo over the volume probed by this survey, 
we expect v v \ = 0.5. We find that eight of our nine fields are within la of 0.5 for sdM/esdM 
with Vtan > 200 km s _1 , and conclude that, as a whole, the survey is complete statistically. 



4.1. Adopted Color-Magnitude Relations 

We use the photometry and parallax data compiled by Gizis (1997) to determine color- 
magnitude relations appropriate for each spectral class of star (M V; sdM; esdM). There are many 
sdM and esdM with accurate distance determinations available, mainly due to the efforts of the 
USNO CCD parallax program ( Monet et al. 1992| ), which allows us to determine the following 



relations: 

M sdM = 4 24 + 2 AQ x — J) (11) 
M esdM = 4 58 + 2.92 x (V — I) (12) 

These relations are valid for (V — I) > 1.6. Our survey is aimed at the fainter, redder M subdwarfs, 
but for comparisons to other surveys of higher mass stars it is useful to have an estimate of the 
main sequence for the bluer subdwarfs. For these stars, the above linear relations are no longer 



appropriate, due to an inflection in the main sequence at V — I ~ 1.5 (c.f., D'Antona & Mazzitelli 



1996| ; |Baraffe et al. 1997| ) . Unfortunately, there are few parallax K subdwarfs and fewer with 



red (RI) photometry. The Hipparcos mission, however, has determined parallaxes for many G 
subdwarfs. We therefore fit a linear relation to the Hipparcos G subdwarfs supplemented by the 
available K subdwarfs from the Fourth Edition of the Yale Parallax Catalog (van Altena et al. 
1996). In this region, models predict that for [m/H] <J —1 the absolute magnitudes are relatively 
insensitive to metallicity ( Baraffe et al. 1997| ). We select only stars with vtan > 220 km s" 1 to 



ensure a halo sample. The data are shown are in Figure §. There are 29 subdwarfs with — < 0.2, 
5 < My < 10, and 0.5 < V - I < 1.4. We find 

M esdG/sdG = 2 9g + 3 lg x _ /) ( 13 ) 

This slope is consistent with the slope measured in globular clusters by Santiago et al. (1996). 

The above relations can only be used for the NGP field, where we have V and I photometry. 
Unfortunately most sdM and esdM lack R photometry, while we have only R and I data for the 
remaining fields in our survey. Hence we have determined the relationship between (V — I)c and 
(r — i)p using Equation |2] and BesselPs (1990) VRIc photometry of nearby stars. We find: 

(r-i) P = 0.1477 + 0.4520(U - I) c - 0.3932(U - lf c + 

0.3606(U - 7)^ - 0.1020(U - I)% + 0.009343(U - I)% (14) 

Applying this transformation to the (V — I)c color of each parallax star, we find that 

Mf dM = 5.37 + 3.53 x(r-i) P (15) 
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M esdM 



5.96 + 4.29 x (r-i) P 



(16) 



Relations 15 and are valid for (r — i)p > 0.77. It is possible that the R-I,V-I transformation 
for the extreme M subdwarfs is not the same as that for disk M dwarfs. While the currently 



available photometry (summarized in Gizis 1997) is unable to resolve this issue, we have obtained 
spectrophotometry ( Gizis 1998a| ) which indicates that there is no significant color term for the 
sdM stars or the earliest (esdK7; esdMO.O) extreme subdwarfs. 



For the G subdwarfs, we transform Equation 13 and obtain: 

M esd/sdG = + g 2 g x , _ y p 



(17) 



Equation 17 applies for 0.31 < (r — i)p < 0.66. 



4.2. Kinematic Corrections 

The luminosity function computed using the techniques described above, although correct in 
shape, represents only the fraction of the halo that has vtan > v cutoff km s . Given a kinematic 
model of the halo, we can calculate the discovery fraction (%) for each field using Monte Carlo 
simulations. BC have shown that variations of ±5 km s _1 in the kinematic parameters ((V),au, 
avi&w) imply variations of less than 4% in sadly, published data for halo stars span a much 
larger range of kinematics. We therefore compare the observed kinematics of our survey stars to 
the predictions of representative halo models. 

Analyses of local (d< 1 kpc) halo stars generally derive mild prograde rotation (e.g., 37 ± 10 
km s , Norris 1986), although recent studies based on more distant (>5 kpc) halo stars have 
found evidence for a retrograde rotating halo (see the review of Majewski 1993). We compare 
our data to the models of Bahcall & Casertano (1986, BC), Norris (1986, N), Casertano et al. 
(1990, CRB), Layden et al. (1996, L), and Beers & Sommer-Larson (1995, BSL). The last is of 
considerable interest as it features a retrograde rotating halo. 

For each of these models, we generated a Monte Carlo simulation of our survey with ~ 5000 
detected stars. The halo velocities are represented by Gaussian distributions in U,V, and W.[] In 
Table ||, we compare the predictions of the U,V, and W velocity distributions of the models to the 
actual data for the 35 metal-poor stars with Mj > 8.5. We emphasize that these are the predicted 
velocity characteristics of the kinematically-biased observed sample, not the true halo kinematics. 
We also give the probability that the observed data and the model simulations are drawn from 
the same distributions, as determined by Kolmogorov-Smirnov tests. The comparison is most 
sensitive to the observed V velocities, which are most important in determining x- R clear that all 



5 We use the standard notation of (U,V,W) for the space velocity components, in which U represents the motion 
towards the Galactic center (I = 0,6 = 0), V represents motion in the direction of Galactic rotation, and W represents 
motion perpendicular to the Galactic plane. 
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the published models agree fairly well with the observed U and W velocity distribution, although 
the N model is significantly poorer for U, due to its smaller velocity dispersion {ajj = 131 km 
s -1 ). The V velocities favor prograde halos, with the BC model (v ro t = 66 km s _1 ) yielding too 
small a mean V velocity and the N model (v ro t = 37 km s— 1) yielding too large a mean velocity. 
The BSL model's retrograde halo is much less likely than the prograde models. Of the published 
models, the Layden et al. model gives the highest probability when all three velocity components 
are considered. 

Guided by the comparison above, we have constructed a composite kinematic model, taking 
the parameters derived by Norris as a starting point and increasing both ajj to match other 
studies, and the halo rotation to a value midway between BC and N. The resulting parameters 
are (v ro t, au, avi&w) = (+50, 140, 106, 85) in km s _1 and we refer to this as our model A. These 
changes are within 1.5a of Norris's estimated uncertainties. The values of 1/x for the BC and N 
models lie within ±7% of the composite values, which are listed in Table |j. Values of 1/x based 
on the Layden et al. halo model are on average 20% smaller. 

We have also calculated models including the IPII (BSL, CRB). Only one in ~ 400 IPII 
stars will enter our vt an > 220 km s -1 sample (although we caution that this number depends 
upon the assumption of Gaussian kinematics). Since the IPII consists of approximately 2 — 5% of 
stars locally, it will contribute less than 10% of our observed halo sample and thus will not bias 
significantly our measurements. This can also been seen in Table 1 of BC, where it is shown that 
the kinematic properties of stars in Eggen's (1979a, 1980) proper motion catalog with vt an > 200, 
vtan > 220 and vtan > 250 km s" 1 are indistinguishable. We expect that the IPII contribution is 
much less than this, since direct spectroscopic observation shows that most IPII M dwarfs both 
locally (nearby stars with 100 < vt an < 220 km s , Gizis 1997) and at height of a few kiloparsecs 
( Rcid et al. 1997| ) are not classified as sdM or esdM. We reject all spectroscopic M V stars from 
our sample, and therefore we expect negligible IPII contamination. 

For each field, we calculate the luminosity function using the Equations [5] to |7] using only 
stars with vt an > 200 km s~ 1 . We combine data for all fields to derive our best estimate luminosity 
function, calculated for 0.5 magnitude bins, which is listed in Table ^. 

These conservative selection criteria ensure a pure Population II halo sample but suffers 
from the uncertainty of the kinematic correction. Since we have obtained spectroscopy for all 
of our reduced proper motion selected stars, we have information on the actual space density of 
metal-poor very-low-mass stars with lower velocities (vtan > 100 km s _1 for esdM; vt an > 125 
km s _1 for sdM). We have therefore computed the space density of sdM and esdM separately, 
using the appropriate velocity cutoff for each metallicity class. The completeness corrections, — , 
for these velocity cutoffs are only ~ 10 — 20%, since virtually all halo stars have vt an > 125 km/s. 
The resulting luminosity functions for the sdM and esdM are listed in Table [7|. This luminosity 
function is remarkably consistent with that estimated from the corrected "pure" halo sample. 
The normalizations differ by no more that 20% from that predicted by our composite model. We 
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discuss these low-vtan samples in more detail in Section |5[ 



5. Discussion 

Our measurement of the Population II luminosity function is aimed at only the limited range 
in luminosity corresponding to the fainter M subdwarfs. Nonetheless, these results offer a means 
of resolving the difference, noted in the Section [j], between Gould et al.'s HST starcount analysis 
at ^ 5 kpc and Dahn et al.'s local trigonometric parallax sample. Before discussing that issue, we 
consider whether our results can be combined with previous measurements of the bright end of the 
halo luminosity function. Since previous studies of the halo luminosity function have presented 
$(My) for all metallicities combined, we have transformed our data using the relation 

$(M V ) = $(M / )^ (18) 



and the color-magnitude relations given in Section 4.1. The results for the combined sdM and 
esdM luminosity function is plotted in Figure [6|. The other luminosity functions shown are the 
local trigonometric parallax samples of Schmidt (1975) and Dahn et al., the local photometric 
parallax sample of Bahcall & Casertano (1986), and the Gould et al. HST starcount analysis. 
We have not plotted Dawson's (1986) statistical analysis of the LHS stars, since it is superceded 
by Dahn et al.'s calculations, which have the advantage of improved data for individual stars. 
Schmidt's results have been updated using literature V-band photometry, but there are only 0-2 
stars per bin, except for the last two bins which have five stars apiece. It should also be noted that, 
except for the HST data, these estimates of $(My) are based on kinematically-selected samples. 
Dahn et al's vt an correction factors, x> agree to within 5% with the values predicted by our model 
A, and we have also used that model to adjust Schmidt's data. The Bahcall & Casertano $(My) 
plotted is derived using their published correction factors. Adopting \ appropriate to model A 
reduces the inferred space densities by 15%, increasing the discrepancy with respect to the other 
studies (including our own results). 

It is obvious from Figure |6| that there are potential problems involved in combining data from 
the several investigations to derive a halo luminosity function spanning the full range in absolute 
magnitude. Of the surveys which include brighter stars near the turnoff, Schmidt's luminosity 
function (as noted originally by Schmidt!) has too few stars to measure reliably the shape of 
the luminosity function. The BC dataset has sufficient stars, but there is a clear, systematic 
offset of almost a factor of two between their derived space densities and Dahn et al's results. 
Bahcall & Casertano's analysis is based on Eggen's (1979a, 1980) survey of southern (5 < 30) 
stars with /j, > 0".7 yr _1 and V< 15, but a comparison of both the relative numbers and V/V max 
distributions of such stars shows no evidence for any significant increase in incompleteness at 
southern declinations. On the other hand, the color-magnitude relation (from Eggen, 1979b) 
which BC use to estimate distances is a poor description of the lower Pop. II main sequence 
(My ^ 9 — 10), and therefore that portion of <&(My) should probably be excluded. The calibration 
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is likely to be more reliable for the more luminous stars, since the main sequence is consistent with 
the Hipparcos data in Figure || for the G subdwarfs. 

One notable aspect of the Eggen/BC sample is that it may be incomplete in the Galactic 
Plane, since it is difficult to identify proper motion stars against the dense stellar background. 
We have simulated halo samples using the BC, A, and L models for both "all-sky" (5 < 30) and 
"non-plane" (5 < 30 and |6| > 20) selection. The kinematic correction factors (1/x for vt an > 220 
km s _1 ) for the "non-plane" samples are greater than the "all-sky" samples by between 8% and 
12% depending on the model. However, our analysis of the Eggen (1979a, 1980) data shows that 
the derived halo space density for stars with \b\ > 20 is actually ~ 40% higher than that derived 
for the entire Eggen sample. This suggests that the BC luminosity function is underestimated by 
~ 30% due to incompleteness in the Galactic Plane. This would account for about a third of the 
discrepancy evident in Figure [6|. 

Naively combining the higher-luminosity stars from BC with Dahn et al. and our own results 
for fainter stars suggests a rather steep luminosity function. Globular clusters provide a potential 
means of testing this result. Figure ^ compares <]?(M/) from the field against results for two of the 
better-studied clusters, using data from Piotto et al. (1997) adjusted to be consistent with the new 
distances derived by Reid (1997) on the basis of main-sequence fitting. That comparison suggests 
that the BC results may underestimate the space densities of stars with 6.5 < Mj < 8, consistent 
with the offset between the BC and Dahn et al results plotted in Figure ||. It should of course 
be noted that the field and cluster luminosity functions may differ - the present-day luminosity 
(mass) function of globular clusters could be affected by dynamical evolution, or alternatively 
the initial mass function of field and cluster stars could be different. Extension of the Dahn et 
al. sample to brighter absolute magnitudes (e.g., expanding the Schmidt sample) should provide 
the most effective means of confirming the shape of the field luminosity function. The Dahn et 
al. luminosity function, transformed via Equation 12 is also plotted, but the transformation is 
uncertain since we have assumed all stars to be esdM, which is not the case — and in any case, 
this preliminary Dahn et al. luminosity function represents a range of metallicities, so it may 
not be fair to compare it to the globular clusters. Nevertheless it is clear that the Dahn et al. 
luminosity function peaks at a fainter absolute magnitude than the globular clusters. 

As noted above, Figure |6| suggests that the Gould et al (1998) HST analysis underestimates 
the local density of late-type halo subdwarfs by a factor of four. To underscore that point, Figure || 
plots our luminosity function against a <1>(M/) representation of the HST data that was computed 
by Gould (private comm.) using the same data and techniques used to derive the My luminosity 
function of Gould et al. (1998). While our study is limited to a restricted range of luminosity, the 
space densities, derived from a completely independent sample of nearby subdwarfs, are clearly 
strongly supportive of the LHS-star analysis. A possible resolution of the discrepancy between 
the distant-halo and local-halo analyses may lie in the composite nature of the galactic halo. 
Sommer-Larsen & Zhen (1990) estimate that at least 40 percent of the local subdwarfs reside in 
a highly-flattened component, which would essentially be absent from the HST dataset. Taken 
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in isolation, this additional component can account only partially for the factor of four offset in 
the derived densities. However, its presence underlines the substantial extrapolation involved 
in transforming the observed number density at large distances from the Plane to a local halo 
subdwarf density. An additional contribution could be due to the choice of local calibrating 
subdwarfs for the Gould et al. analysis, which may be biased towards higher metallicities. This 
would lead to overestimated My and distance, and hence an underestimate of the space density. 
For a discussion of possible systematic errors, see Gould et al.'s discussion. 

The kinematics of the very-low-mass halo stars are consistent with those found for higher 
mass metal-poor stars. The strongest conclusion is that prograde rotating halos are favored. If the 
expanded halo sample that includes lower tangential velocity stars is included, then kinematics 
like those of model A is favored over Layden et al.'s RR Lyrae kinematics, since the expanded 
luminosity function is more consistent with the model A corrected pure halo sample (Tables || 
and 0). There may be an excess of sdM at low velocities over the halo model A predictions. Since 
metallicities of [m/H] ~ —1.2 ± 0.3 are expected to have some IPII contribution, this is not too 
surprising. If the two best determined bins are used, then the excess is only 19%. Chiba & Yoshii 
(1998) find an excess of ~ 20% in the range —1.4 < [Fe/H] < —1.0 based upon Hipparcos data 
for red giant and RR Lyrae stars, in excellent agreement with our results. 

The mass function can be estimated using the observed luminosity function and a 
mass- luminosity relation as follows: 



Since there are no direct determinations of mass for any sdM or esdM stars, we must rely upon 
theoretical mass-luminosity relations. These must be viewed with caution - for example, Reid & 
Gizis (1997) have shown that, in the case of the Baraffe et al solar-abundance models, the slope 
deduced from nearby-star data for the disk mass function depends on the passband used in the 
analysis. Our data do not provide sufficient range in Mj to estimate the mass function reliably, 
even given a good mass- luminosity relationship, but it does allow us to estimate the relative 
numbers of sdM and esdM stars at a mass of 0.2M Q . We transform the observed pure halo Mj 
luminosity functions into mass functions using the Baraffe et al. (1997) models. We adopt the 
[m/H] = —1.0 and [m/H] = —1.5 interior models for the sdM and esdM respectively, and fit a 
power-law to the mass function in order to obtain the best estimate of the numbers. The ratio 
of sdM to esdM at 0.225M Q is 0.38. In comparison, we obtain a ratio of 0.36 if we adopt the 
Alexander et al. (1997) models, and 0.60 for the D'Antona & Mazzitelli (1996) tracks — but the 
latter match the observed V — I, Mj main sequence poorly. These results are consistent with the 
Carney et al. (1994) distribution of metallicities for higher mass stars (mainly G subdwarfs) - 
41% of their stars in a pure (V < —220 km s _1 ) halo sample have [m/H] > —1.5. Thus, the 




(19) 



^f(M) = oc M 

V ' dM 



—a 



(20) 



-16- 



metallicity distribution for the very-low-mass stars of the halo (M ~ 0.2M©) is apparently the 
same as that of stars on the upper main sequence (M ~ 0.7M Q ). Barring a cosmic conspiracy 
of offsetting parameters, the most reasonable explanation is that the mass function in the halo 
is independent of metallicity. Although our data is inadequate to constrain the mass function 
well, we note that the best-fit parameters are a ~ 0.5 ± 1.6 for the esdM and 1.3 ± 0.9 for the 
sdM. These compare with a power-law index of a = 1.05 ± 0.15 for disk stars in the mass range 
1.0 > > 0.1 ( Reid fe Gizis 19971 ). 

The consistency of the metallicity distribution and kinematics of the low-mass stars (sdM 
and esdM) is not surprising, but nevertheless important. Studies of the halo, whether using 
evolved stars such as RR Lyrae or red giants, or main sequence stars (F,G subdwarfs), assume 
that the properties of stars in the narrow mass range from just below the turnoff to the tip of the 
asymptotic giant branch (0.6 — O.8M0, are representative of the halo population as a whole. Our 
results verify one aspect of this assumption for stars near the bottom of the halo main sequence. 

6. Summary 

We have derived a new estimate of the space density of very low mass metal-poor stars. 
Comparisons with other surveys shows that we are in agreement with the Dahn et al. (1995) 
survey based upon trigonometric parallaxes of LHS catalog stars. This implies that the local space 
density of metal-poor stars is higher than that predicted by HST observations of more distant 
stars. 

Our work provides measurements of the kinematics and metallicity distribution of very-low- 
mass stars. These measurements are consistent with those from studies of higher-mass stars in the 
halo. We find that the esdM are 2.5 times as common as sdM at ~ 0.2M Q . This is consistent with 
the relative numbers of metal-poor ([m/H] < —1.5) and metal-rich ([m.H] > —1.5) G subdwarfs 
found by Carney et al. (1994). Note that the two metallicity scales are comparable because Gizis 
& Reid (1997) have found agreement between the Carney et al. metallicities of FGK subdwarf 
primaries with metallicities based upon G97. The agreement of the metallicity distribution at 
both low and high mass implies that the mass function is not a strong function of metallicity for 
(M ^ O.7M0) halo stars. This result is supported by the general similarity of the mass functions 
derived for the two metallicity bins considered; however, those mass functions depend upon 
uncertain model mass- luminosity transformations. 
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of Palomar Observatory for their capable support. We also thank M. Schmidt for useful 
conversations and A. Gould for providing the I-band HST luminosity function. J. E.G. is 
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A. Stars of Particular Interest 

Most of our subdwarfs lie at distances of greater than 200 parsecs, and therefore are more 
difficult to study than LHS Catalog subdwarfs. Three of our stars, however, are close enough that 
additional observations may be profitable. All were previously identified by Luyten (1979-1980). 
The positions (B1950, at the epoch given in Table p]) and other data for these stars are given in 
Table ||. Spectra of the subdwarfs LP 622-7 and LP 382-40 are plotted in Figure ||. LP 622-7 has a 
proper motion too large to be included in our complete sample, and was observed at the Palomar 
60-in. in June 1995. The radial velocity may be unreliable due to a problem with the calibrating 
arc lamp exposure. 

LP 589-7 is one of the coolest subdwarfs known. Its spectral indices are Ti05= 0.68, 
CaHl=0.48, CaH2=0.27, and CaH3= 0.50. In Figure [l0|, we show that this star is extremely 
similar to the esdM5.0 QGizis 1997Q star LHS 3061 (Mj = 11.65, |Monet et al. 1992|) . The 



conspicuous TiO absorption at 70504 distinguishes both these stars and the slightly cooler 
LHS 1742a (Mj = 11.69, esdM5.5) from the other two very cool extreme subdwarfs LHS 205a 



{Mi = 11.65, esdM5.0) and LHS 1826 (Mj unknown, esdM6.0; |Gizis fc Reid 1997|) . LP 589-7 is 



just as close as LHS 1742a, so it is an excellent candidate for a trigonometric parallax. 

There is a hint of Ha emission in our spectrum of LP 589-7. The feature, which may not real 
but lies at the correct wavelength, has an equivalent width of ~ 1.2 A. Two sdM systems (LHS 
482 and LHS 2497) are known to show emission due to close companions flGizis 1998b ). Since LP 



589-7 is near the hydrogen burning limit, any companion must be either a very faint, old white 
dwarf (as for the LHS 482 system), a nearly equal luminosity M companion (as for LHS 2497), 
or else a brown dwarf (if it proves to be too faint to contribute detectable light). Alternatively, 
LP 589-7 could a "young" metal-poor star with an age of up to a few billion years. If an SB2, 
its estimated distance and tangential velocity should be increased by up to 40%. All this is 
speculative, however, since the Ha feature is weak and may not be real. 
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Table 1. Survey Fields 
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Table 2. Relative to Absolute Proper Motions 
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Table 3. Spectroscopic Indices 
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Table 4. Kinematic Corrections 
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17.3 


2.30 


NGP 


28.00 


17.0 


2.01 


868 


25.29 


17.3 


2.11 


513 


25.29 


17.3 


2.18 


889 


25.30 


17.3 


3.70 


890 


25.26 


17.3 


3.56 



-23- 



Table 5. Kinematic Comparison 





Data 


BC 


N 


L 


CRB 


BSL 


A 




171 


171 


146 


179 


165 


155 


158 


Probability (U) 




0.43 


0.15 


0.55 


0.29 


0.19 


0.31 


(V) 


-234 


-222 


-244 


-240 


-252 


-264 


-233 


ay 


76 


82 


82 


84 


73 


76 


85 


Probability (V) 




0.37 


0.49 


0.75 


0.22 


0.05 


0.87 




90 


77 


86 


97 


94 


105 


86 


Probability (W) 




0.38 


0.49 


0.32 


0.38 


0.28 


0.49 



Table 6. Pure Halo Luminosity Functions 





$ [A] a$ [A] $ [L] 


o-$ [L] 


N 




10~ 5 pc~ 3 Mag" 1 10~ 5 pc 


~ 3 Mag" 1 






sdM (vtan > 200 km s' 






8.5 


2.2 0.7 1.8 


0.5 


16 


9.0 


2.0 0.6 1.6 


0.5 


11 


9.5 


2.0 0.9 1.6 


0.8 


5 


10.0 


2.0 2.0 1.6 


1.6 


1 




esdM (vtan > 200 km s 


- 1 ) 




9.5 


4.6 1.8 3.9 


1.5 


10 


10.0 


3.6 1.6 2.9 


1.3 


5 


10.5 


3.1 2.3 2.6 


2.0 


2 
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Table 7. Expanded Halo Luminosity Functions 



Mi $ [A] a$ [A] $ [L] g <j> [L] N 
10~ 5 pc~ 3 Mag" 1 1(T 5 pc- 3 Mag" 1 



sdM (vtan > 125 km s" 1 ) 



8.5 


2.6 


0.8 


2.5 


0.7 


20 


9.0 


2.4 


0.8 


2.2 


0.7 


14 


9.5 


4.4 


1.3 


4.1 


1.2 


12 


10.0 


3.1 


1.6 


2.9 


1.5 


4 


10.5 


3.6 


3.6 


3.4 


3.4 


1 




esdM (vtan 


> 100 km s 


- 1 ) 




9.5 


6.9 


1.8 


6.6 


1.7 


20 


10.0 


2.5 


1.0 


2.4 


1.0 


6 


10.5 


4.1 


2.1 


3.9 


2.0 


4 


11.0 













11.5 













12.0 


15. 


15. 


15. 


15. 


1 



Table 8. Nearby Subdwarfs 



Name 


a 


S 


Ma 


AM 


i P 


r — i p 


M, 


d 


Vtan 


Vrad 


Sp. Type 


LP 589-7 


28.72333 


+1.03647 


-0.12 


-0.33 


15.90 


1.38 


11.77 


64 


106 


-65 


esdM5 


LP 382-40 


223.53375 


+23.83433 


-0.28 


-0.09 


16.47 


1.08 


10.59 


150 


210 


-230 


esdM2.5/3 


LP 622-7 


227.18958 


-1.02281 


+0.20 


-0.43 


13.87 


0.87 


9.69 


69 


155 


300 


esdK7/M0 
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Fig. 1. — A color - reduced proper motion (V — I, Hy) diagram for a nearly volume limited sample 
of stars within 25 parsecs. The 2111 observable known nearby M dwarfs (Reid et al. 1995; Hawley 
et al. 1996) are plotted. For some stars, the V-I photometry was estimated on the basis of their 
spectral type. The solid line plots a disk population using the Stobie et al. (1989) main sequence 
and vtan = 50 km s" 1 . The dotted line plots the same main sequence with vtan = 220 km s -1 . The 
dashed line plots a main sequence that is two magnitudes subluminous, roughly corresponding to 




Fig. 2. — Comparison of magnitudes for matches within 3 arcsec for the POSSII field 513 and 
POSSI plate E102. The outlier points are due to galaxies and merged sources which are treated 
differently by the two measuring machines. 
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Fig. 4. — The reduced proper motion diagram for Field 513. Halo stars lie below the dotted line, 
which corresponds to 220 km s _1 for near-solar metallicity disk stars. 
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Fig. 5. — The HR Diagram for stars using trigonometric parallaxes. Crosses are stars earlier than 
type M with vtan > 220 km s _1 from the Yale Parallax Catalog ( |van Altena et al. 1995| ) or the 
Hipparcos Catalog, solid squares are sdM (Chapter 2), and solid triangles are esdM. The disk 
sequence is shown using nearby stars (open squares) 
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Fig. 6. — Comparison of our halo luminosity function with previously published halo V-band luminosity 
functions. Our luminosity functions, measured in Mj have been transformed to $(Afy) using Equations [l2| 
(esdM) and [ll] (sdM) and then added together (solid squares). The luminosity functions shown are Schmidt 
(1975, open squares with dashed line), Bahcall & Casertano (1986, solid circles with long-dashed line), Dahn 
et al. (1995, open squares with solid line), and Gould et al. (1998, asterisks). Note that our luminosity 
function supports the higher value measured for local LHS stars by Dahn et al. rather than the lower value 
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Fig. 7. — The globular clusters NGC 6341 (open triangles) and NGC 7099 (open circles) compared to the 
csdM luminosity function (solid triangles), Dahn et al. (open squares), and Bahcall & Casertano (solid 
circles). The BC luminosity functions has been transformed to $(Mj), no account has been taken of any 
metallicity spread; BC's original kinematic correction has been applied. Correcting for the latter two effects 
would shift the BC luminosity function downwards. When the globular clusters are normalized to our 
luminosity function, their luminosity is parallel to but above the BC luminosity function for brighter stars. 
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Fig. 8. — Our estimated luminosity function for sdM (solid triangles), esdM (solid squares), and 
both combined (open squares). Also shown is the Gould et al. HST luminosity function (asterisks). 
Due to crowding, the error bars have been suppressed for the individual esdM and sdM datasets. It 
is clear that the HST luminosity function predicts substantially fewer stars than are actually seen. 
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Fig. 9. — The nearby subdwarfs LP 622-7 and LP 382-40. The apparent emission lines at the blue 
end of the spectra are due to cosmic rays. 
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Fig. 10. — The nearby subdwarf LP 589-7 is nearly identical to LHS 3061. Both are classified as 
esdM5.0 using the Gizis (1997) system. Note the TiO absorption in both stars at 7050A, which 
distinguishes them from the other known esdM5.0 star, LHS 205a. 



